The main code of the work is in 'CSTR.py', 
including the CSTR model, the optimization, and the plot. 
The code can be modified with different parameter settings, 
including the initial condition, the tuning parameters of the shrunk target zone, etc.
The computation code of the CIS can be found in the folder "cis_song_toMia", with credits to Ben. 
The CIS boundary info files can be found in the folder as well.